Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 3 March 2013 (MN M^jX style file v2.2) 



Diagnostics of Baryonic Cooling in Lensing Galaxies 



(N 

o 

< 

in 
(N 

O 

u 

43 
6 



> 

in 
o 

(N 



IS 



Dominik Leier 1 *, Ignacio Ferreras 2 , Prasenjit Saha 3 

1 Astronomisches Rechen-Institut, Zentrum fur Astronomie der Universitdt Heidelberg, Monchhofstrasse 12-14, 69120 Heidelberg, Germany 

2 Mullard Space Science Laboratory, University College London, Holmbury St Mary, Dorking, Surrey RH5 6NT, UK 

3 Institute for Theoretical Physics, University of Zurich, Winterthurerstrasse 190, CH-8057 Zurich, Switzerland 



Accepted 2012 April 24. Received 2012 February 27; in original form 2011 September 26. 



ABSTRACT 

Theoretical studies of structure formation find an inverse proportionality between the 
concentration of dark matter haloes and virial mass. This trend has been recently con- 
firmed for M vir > 6 x f 12 Mq by the observation of the X-ray emission from the hot 
halo gas. We present an alternative approach to this problem, exploring the concen- 
tration of dark matter haloes over galaxy scales on a sample of 18 early-type systems. 
Our c — M vir relation is consistent with the X-ray analysis, extending towards lower 
virial masses, covering the range from 4 x IO^Mq up to 5 x 10 12 M q . A combination 
of the lensing analysis along with photometric data allows us to constrain the baryon 
fraction within a few effective radii, which is compared with prescriptions for adia- 
batic contraction (AC) of the dark matter haloes. We find that the standard methods 
for AC are strongly disfavored, requiring additional mechanisms - such as mass loss 
during the contraction process - to play a role during the phases following the collapse 
of the haloes. 

Key words: gravitational lensing - galaxies: elliptical and lenticular, cD - galaxies: 
evolution - galaxies: haloes - galaxies: stellar content - dark matter. 



1 INTRODUCTION 

Dark matter haloes constitute the scaffolding on which the 
luminous component of cosmic structure can be detected in 
the form of galaxies. The connection between ordinary mat- 
ter (i.e. "baryons" ) and the dominant dark matter can only 
be found indirectly because of the elusive nature of the lat- 
ter. Theoretical studies give useful insight under a number of 
assumptions about the properties of the dark matter compo- 
nent. N-body simula tions predict a universal density profile 
|Navarro et al.l ll996) driven by two parameters, the mass of 
the halo - usually defined out to a virial radius - and the 
concentration, given as the ratio between the scale length of 
the halo (r s , where Alogp/Alogr = —2) an d the virial ra- 
dius. A number of N-body simulations (see e.g. lBullock et afl 
l200ll : lNeto et al.ll2007l ; iMaccio et alj|2007h display a signifi- 
cant trend between these two parameters, mainly driven by 
the hierarchical buildup of structure. Massive haloes assem- 
ble at later epochs - when the background density is lower 
because of the expansion of the Universe. Hence, we ex- 
pect a trend whereby concentration decreases with increas- 
ing galaxy mass. An observational confirmation of this the- 
oretical result by means of model-independent mass recon- 
struction gives insights in the interplay between luminous 
and dark matter. 
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Over cluster scales it is possible to explore the dark 
matter halo via its effect on the gravitational potential. The 
X-ray bremsstrahlung emission from the hot gas in the in- 
tracluster medium acts as a tracer of the potential. If as- 
sumptions are made about the dynamical state of the clus- 
ter, o n e can constrain the halo properties (see e.g. lSato et al.l 
2000). lBuote et all (|2007h found a significant correlation be- 
tween dark matter halo and mass, as expected from theo- 
retical studies. Their sample covers a range from massive 
early-type galaxies up to galaxy clusters. In this paper, we 
extend the observational effort towards lower masses, con- 
straining the haloes over galaxy scales by the use of strong 
gravitational lensing. By targeting a sample of strong lenses 
at moderate redshift (z~0.5) we probe the mass distribu- 
tion out to a few (~4) effective radii. Other approaches to 
probe dark matter haloes over galaxy scales involve the use 
of dy namical tracers such as the bulk of the stellar popula- 
tions jGerhard et alj|200lh . While Sauron data, as used in 
(|2006l ) is restricted to 1-2 effective radii 
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by the surface brightness detection limit of the observa- 
tions and contains thus little information on the proper- 
ties of the halos, studies based on mo re extended data are 
available (see e.g. iThomas et al.l 12009). Addition ally, plan- 
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matter profile. Being evolved phases of the underlying stellar 
populations, planetary nebulae can be consi dered unbiased 
trace rs of the gravitational potential (see e.g. ICoccato et al.l 
2009). Through their emission lines, it is possible to trace 
their kinematics out to l arge distances, reaching out to 
~ 5fl. |Douglas et all 120021 ). However, the interpretation of 
the results is difficult because of the uncertainties regarding 
the parameteri sation of the halo mas s, anisotropy, shape, 
or inclination (Ide Lorenzi et all I2009T ). Gravitational lens- 
ing studies do not suffer from the inherent degeneracies of 
methods regarding the modelling of the dynamical tracers, 
although it is fair to say that lensing studies have other mod- 
elling degeneracies, as we discuss in the following Section. 
Ultimately, a comparison between all these methods is key 
to a robust assessment of the c — M vir relation. 

Our recent study of stellar and total mass in lensing 
galaxies l|Leier et aUl20lJ . hereafter LFSF) indicated an in- 
verse trend of concentration with mass. Those results, how- 
ever, applied to radii much smaller than the virial radius r vir . 
In this work, we extend this analysis by extrapolating the in- 
ferred dark-matter profiles out to r vir , to determine whether 
the concentration/mass trend persists to lower masses, i.e., 
galaxy haloes. We then try to reconstruct the possible con- 
centrations of the dark-matter haloes before adiabatic con- 
traction due to the baryons. 

We consider a sample of 18 early-type lensing galax- 
ies. In LFSF these galaxies (plus two disk galaxies and one 
ongoing major merger) were decomposed into stellar and 
dark-matter profiles. The stellar-mass or M ste i profiles were 
obtained by fitting stellar population synthesis models to the 
star light. The lensing- mass or Mi ona profiles were obtained 
from lens models. The difference between these is assumed 
to be the dark-matter profile. 

In Section [2] of the present paper we summarize the 
method of deriving the dark- matter profile as Mi 8ns — M at ei- 
We also consider the technique of simply fitting a two- 
component lens model. The latter technique, in a test case 
(see Fig. [XJ) appears adequate for estimating Mi eM but sig- 
nificantly overestimates M Bte i. 

In Section[3]we fit well-known NFW and Hernquist pro- 
files to the dark and stellar-mass profiles respectively. Fig- 
ure [2] shows the NFW and Hernquist parameter estimates 
and uncertainties for two of the galaxies, while Figure 
shows the dark-matter profiles for the same two galaxies. 
The NFW fits automatically provide a virial mass M vir and 
a concentration c, in effect extrapolating the dark-matter 
profile out to the virial radius r vir . Figure [4] shows M vir 
and c for all 18 galaxies. The trend shown in iBuote et al.l 
l|2007l ) is seen to extend down to virial masses of 10 12 M Q . 
We remark that the NFW fits shows a characteristic banana- 
shaped near-degeneracy between the virial mass M vir and 
the concentration c. These contribute a spurious inverse cor- 
relation between M vir and c, but they are much smaller than 
the overall trend. 

In Section [4 ] we use abun dance matching (e.g., 
iMoster et al"1l201oTlGuo et alj|20ld ) to derive a virial mass 
M* r M directly from M atol . The two estimates M* r M and M vir 
tend to agree in the majority, but there are cases of strong 
disagreement. Interestingly, the latter are all galaxies in 
dense environments. We also consider the option of con- 
straining the NFW fit such that M vir = A/* r M . Figure [5] 
shows how the mass profiles get modified if this is done, 



while Figure [6] shows how the c-M vir distribution changes. 
In the latter case, the scatter increases considerably. 

In Section [S] we attempt to reconstruct the initial con- 
centrations, by fitti ng the adiabatic-contraction model of 
iGnedin et ail (|2004 ). We find that the usual prescriptions 
for adiabatic-contraction imply unrealistically low values of 
c init , but weaker adiabatic contractions do fit our results (see 
Figure [7]). By tweaking the average radius in the adiabatic- 
contraction prescription (which can be interpreted as mass 
loss during adiabatic contraction) we can obtain agreement 
with the data. The inferred c init are shown in Figure [9l from 
which it appears that adiabatic contraction can explain part 
of the c-A/ vir trend but is unlikely to be the sole origin of it. 



2 MULTI-COMPONENT FITTING VS 
STELLAR POPULATION MODELLING 

The starting point of the present work is the models in LFSF 
of the projected stellar and total surface mass density from 
a sample of lensing galaxies. We obtained independent maps 
for the stellar mass and total mass, using archival data from 
the CASTLeS survejQ- The maps of stellar mass, M ato i, were 
derived by fitting stellar population synthes is models t o pho - 
tometry in two or more bands assuming a IChabrierl (|2003l ) 
initial mass function (IMF). The total or lens mass, Mi Gna , 
was mapped by computing pixellated lens models that fit- 
ted the lensed images and (where available) time delays. 
Detailed error estimates were derived in both cases. The en- 
closed total mass is well constrained at projected radii where 
images are present. At smaller and larger radii, Mi ona be- 
comes progressively more uncertain. The outer radius of the 
mass maps is 2_Ri ons where i?i ona is the radius of the outer- 
most image. Since iii cnB depends on the redshift and details 
of the source position, 2i?i ona varies greatly among galaxies 
— from a quarter of the half-light radius (R e ) to several R c . 

Of the sample modelled in LFSF, 18 galaxies are early 
type systems. We exclude the Einstein Cross Q2237, which 
is the bulge of a spiral galaxy; B1600, which is likely to be 
a late-type galaxy viewed edge-on; and B1608, which is an 
ongoing merger. For these 18 early type galaxies, there is 
no evidence of a significant gaseous component, and hence 
we may assume that AM — M\ eas — M atc i is a map of the 
dark matter distribution. The lensing maps tend to have 
similar orientation to the stellar mass, and hence the AM 
maps are fairly elliptical as well (F erreras et al.l l2008t ) . In 
LFSF we obtained enclosed dark matter profiles AA/(< R) 
using a circularized aperture along the elliptical isophotes, 
i.e. following the luminous distribution. 

An alternat ive approach (see e.g. lAuger et al.l l2010l : 
iTrott et al.ll201ol ) consists of fitting a parametric lens model 
with separate components for stellar and dark matter. If 
the stellar component can be correctly recovered by this 
method, the analysis based on stellar population synthe- 
sis constrained by multiband photometry will be dispens- 
able. The two approaches are contrasted in Fig. [T] in the 
case of the quad PG1115+080. On the one hand we pre- 
pared separate models for Mi ona and M ato i — a pixellated 
lens model for Mi cna and a stellar-population model from 

http://www.cfa.harvard.edu/castles 
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Figure 1. Comparison of two different modeling strategies on the lens PG1115+080. Left: A fit to a three-component lens model (stellar, 
dark matter, and external group) without population synthesis. Red dots mark the image positions and the red curve is the model critical 
curve. The blue dot is the model source position and the green curves show the model caustics. The grayscale indicates stellar-mass 
fraction, while the black ellipse indicates the ellipticity and position angle of the stellar component. The semi-major axis of the latter is 
arbitrary, and set here to 2iJ lcnB . Right: The red and cyan dots show the stellar and total enclosed mass respectively from the model in 
the left panel. The red error bars enclosing the red dots correspond to a 1<t region around the best x 2 - The gray and black dots with error 
bars are the stellar mass and total enclosed mass respectively from the models in LFSF, which use pixellated lens models and population 
synthesis. 



the photometry for M ste i. On the other hand, we fitted the 
lensing data to a multi-component parametric lens model: 
a de Vaucouleurs profile, plus an NFW halo, together with 
a singular isothermal sphere, adding external shear to ac- 
count for a nearby ga laxy group. We used the gravlens 
program (|Keetonll200lh together with Markov-chain Monte- 
Carlo (MCMC) to search for the best fit parameters. The ef- 
fective radius was constrained to lie with in the observational 
uncer tainty R c = 0.85 ± 0.07 arcsec (|Treu fc Koopmansl 
2002). The positions, ellipticities and position angle were 
also allowed to vary. 



We see from Fig. [TJ that the parametric and pixellated 
lens models give similar results for the total-mass profile. 
The pixellated method provides uncertainty estimates be- 
cause it generates an ensemble of models. It is also compu- 
tationally faster. However, the stellar mass is strongly over- 
estimated by the parametric lens model, compared to the 
estimates based on population synthesis. In other words, 
the attempt to infer stellar masses from the lensing data 
alone fails. An explanation for this is suggested by the error 
bars on the total-mass. We see that lensing provides a good 
estimate of the enclosed mass at radii comparable to the 
images, but gets progressively more uncertain as we move 
inwards or outwards. Thus, attempting to extract the profile 
of a sub-component from this already-uncertain total-mass 
profile (without adding more data) will tend to amplify the 



Hence, for the rest of this paper, we will use the separate 
models of M lBnB and M Bta i from LFSF. 



3 VIRIAL MASS AND CONCENTRATION 

With the maps of AM = Mi BnB — M 5to i surface mass density 
in hand, we now proceed to estimate a virial mass M vir and a 
concentration c for each of the lensing galaxies. The method 
we adopt is to fit the profiles of AM(< R) to the cumulative 
projected mass of an NFW profile, which is given by: 



'(< R) = 4Tvp s ri x T(R,r s 



M 
where 



,F(r,r s ) = In — + < 

AT a 



(1) 



■ cosh 1 ?f (r < r a ) 



(r = r s ) (2) 



(r > r s 



Here r 3 and p 3 are the scale radius and scale density pa- 
rameters on which the NFW profile depends. For AM(< R) 

a\ u where 5 MlenB is 



we assumed an error a 



half of the 90% confidence interval given by the ensemble 
of lens-mass models, and <XM stel is the standard deviation of 
stellar mass from population synthesis. The best fit values 
of the parameters can be found in Table [2] 

For two example lenses, Q0047-280 and HE2149-274, 
we illustrate the results in more detail in Figs. [2] and [3] The 
upper panels of Fig. [2] show parameter fits and x 2 contours. 
Note that the axes of these two panels are not simply r s 
and p s but rather r s /R e and p s r 2 , in units of MqR~ x . This 
choice tends to illustrate the parameters better. To check 
for the possibility of multiple local x 2 -minima we generated 
MCMC chains with 10 steps for each lens. Such additional 
minima can be excluded for physically interesting parameter 
values. As a further check, we compute as before the NFW 
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Figure 2. Parameter fits to dark matter and stellar (AM and 
M sto i) profiles for the lens galaxies B0047 and HE2149. Top: A 
X 2 map of the NFW parameter space for AM, in grayscale with 
contours of A% 2 = 1,2,3. Red crosses mark the overall best fit. 
Cyan crosses mark best fits to profiles at the steep and shallow 
ends of the confidence region. Bottom: The same for the Hernquist 
parameters for the luminosity L. 



parameters based on a x search that best fit the steep- 
est and shallowest profiles A/i cn3 (< R) allowed by the LFSF 
analysis of the lensing data. These are indicated by cyan 
crosses in Fig. [2] and, as expected, are roughly in the region 
of the fits to AA/(< R). 

The lower panels in Fig. [5] show parameter fits to the 
distributions of luminous matter, for the same two example 
galaxies, that we will use lat e r, in S ection[5l We use the well- 
known profile of iHernoui st (1990). The enclosed projected 
form of the Hernquist (analogous to Eq.[T]for the NFW) is: 



L(< R) 
where 
X(r,r h ) = 



M 
— 



R_\ 2 X(R,r h )-l 
r h ) 



h =, , 2 sech 1 (r/rh) r ^r h 
77 / „ ; sec" 1 (r/r h ) r > r h 



(3) 



(4) 



Again, the parameters and lcr errors as well as our values 
for r- vir are given in Table [2] Note in the figure we show the 
contours with respect to total luminosity, i.e. L = M/T. 

Figure [3] shows profiles of AM(< R), together with 
NFW fits and uncertainties. There is a tendency for the 
innermost point (in these two examples as well in other 
lenses in our study) to be higher than the fit. We note 
that various simulations (|Moore et afl 1 19981 : iNavarro et al.l 
120041 : iDiemand et alj 120051 ) indicate a somewhat steeper 
slope than the original NFW. Recentlv lCardone et al.l (|201 lh 
advocated a generalized NFW profile with an additional pa- 
rameter. Furthermore, the presence of baryons will tend to 
steepen the central dark-matter profile through adiabatic 
contraction (although we note that feedback effects, such as 
baryon ejecta from supernovae-driven winds, or dynamical 



Figure 3. Inferred dark-matter profile AM and NFW fits, for the 
lens galaxies B0047 and HE2149. The red (cyan) line corresponds 
to the red (cyan) crosses in Figure [2] 



interactions with smaller structures could have the opposite 
effect, making the inner dark matter profile shallower). We 
will address the issue of adiabatic contraction later, in Sec- 
tion[4] For now we assume that a projected NFW depending 
on scale radius r s and the normalization p s sufficiently de- 
scribe the data. 

^From the NFW parameters r s , p s and given the red- 
shift of the halo, the virial mass M vir and concentration c are 
easily derived. Consider first the mass enclosed in a sphere 
(not to be confused with the cylindrical enclosed mass Eq. 



M ap h(< r) = 4np B r s < In 



1 + - S )-TTT 



and the mean enclosed density within a given radius 



(5) 



(6) 



The virial radius is the r at which the mean enclosed density 
equals a certain multiple A c of the critical density, namely: 



(P(< fvix)) = A cPc (z) 

and the mass within the virial radius 

4 ., 

M vir = A c p c (z) X -7T< ir 

is the virial mass. The concentration is defined as 

fvir 



c = 



(7) 



(8) 



(9) 



The value for the overdensity is 
A c = 18tt 2 + 82s - 39x 2 
where x = (Q M (1 + z) 3 /E(z 

^13 



(10) 

1 and E{zf = Sl M (l + 



z)~ + S} A l|Brvan fc Normanlll998h . This value of A c gives the 
exact virial radiu s for a top-ha t perturbation that has just 
virialized (see e.g.[Pecblcs 1980]). The galaxies in our sample 
would have virialized well before the observed epoch. Hence, 
if the observed redshift is used to derive an r vlr , the value is 
unlikely to have the dynamical interpretation of a virial ra- 
dius. Neverth eless, since such a definition of ry ir is comm only 
adopted (e.g. iBrvan fc Normanlll998l : iBuote et al1l2007l ) we 
adopt it in the present work. 

The values of c and M vir , along with errors calculated 
according to the projected lcr regions of Figure [2] are quan- 
tities are listed in Table. [2] Figure [4] plots the values — there 



© 0000 RAS, MNRAS 000, 000-000 



Diagnostics of Baryonic Cooling in Lensing Galaxies 5 




Figure 4. Concentration versus virial mass. The red dots rep- 
resent data from this study. The grey circles shows data from 
iBuote et al.l j2007h ■ The grey contours show the la error region 
from Fig. [l] The red dashed line shows a bootstrapping fit to 
our data. The solid red line shows the result of a bivariate fitting 
method for correlated err ors and instrinsic scatter (BCES) by 
lAkritas fc Bershadvl lll996|'l applied to the data. The same method 
was used bv lBuofaTetalJ 1120071) to obtain the dashed grey line. 



have been no previous data showing M vlr vs c to such small 
virial mas ses. For comparison , Figure [4] also shows the re- 
sults from IBuote et al.l ((2007) of the X-ray c -M viT relation 
for 3 9 systems. These ra nge in M vir from 6 x 10 12 to 2 x 10 15 
M e . IBuote et all (|2007h fit a power-law 



Cl4 / M vir 



1 + 2 KMl 



(11) 



where M M = lO 14 h _1 M is a reference mass and C14 and a 
are constants independent of M, obtaining a — — 0.17±0.03. 
Leaving out the (1 + .z)-term gives a — —0.20 ± 0.03. 
As in their analysis we use a bivariate fitting method 
for correlated errors and intrinsic scatter (BCES) due to 
lAkritas fc Bershadvl (|l996h . which gives a = -0.40 ± 0.06. 

Note how in Figure |4]the projected la contours cover a 
smaller region than simple vertical and horizontal error bars 
would imply. This suggests that more information might be 
available with the appropriate tools. Hence, we tried alter- 
native fitting approaches. In a parametric boostrap we ran- 
domly resample over arbitrary points within the la con- 
tours, so that one point per lens is used for an ordinary 
least square fit. This is done for 10 4 realizations. In a later 
run we ease the requirement of one point per lens and pick 
instead a number n < 18 out of the total number of lenses 
to perform the resampling. In all cases the mean value of the 
slope a stays the same, but its distribution gets broader for 
smaller values of n. For n = 18 the bootstrapping analysis 
yields a — —0.42 ±0.08. We also employed a piecewise anal- 
ysis to check how the slope a of the relation evolves and to 
see whether fits in common mass range yield similar results. 
Furthermore we fit a combined sample of 57 objects. The 
results are shown in Table [1] 

Going from high to low M vir the slope increases from 
-0.10 ± 0.05 f or M vir > W 14 M a to -0.20 ± 0.13 for 
M vir < 10 14 M Q (|Buote et al.ll2007h and finally -0.40 ±0.06 



Figure 6. Similar to Fig. [4] but showing how the fits change when 
abundance matching is included. The arrows go from the original 
position (red dots in Fig. [4} to the position f or which the parame - 
ters include information on virial mass from lMoster et al.1 j201(t) . 
The red arrows mark the 10 lenses for which the extrapolated 
NFW analysis and abundance matching give consistent answers. 



for M vil < 4 x 10 13 . For the mass range between 6 x 10 
Mq and 1 x 10 14 Mq, where the two samples overlap, 
we do not find significant differences — a two-dimensional 
Kolmogorov-Smirnov test for the overlapping region gives 
a p-value of ~ 0.5 under the null hypothesis that both 
samples are drawn from the same population. However, it 
should be noted that the reduced sample size and consider- 
able scatter leads to large errors for both sa mples. A trend 



of a w ith virial mass w as first sugges t ed by iNavarro et al 



|l996h and confirmed bv lBullock et al] (|200lh and lEke et al 
(2001) for simulations. Higher normalization factors com- 
pared to simulations are also kn own from a lensing study by 
IComerford fc Nataraianl j|2007h . 



4 COMPARISON WITH ABUNDANCE 
MATCHING 

Thus far we have computed c and M vir under the assump- 
tions that (a) an NFW profile is a good representation of 
the dark matter profile beyond the radial range probed for 
lens galaxies in LFSF, and (b) the dark matter profile is well 
constrained by pixellated studies of stellar and total mass, 
meaning also that the probed radial range is sensitive to the 
scale radius of the NFW profile and that the uncertainties 
give a robust estimate of the suitable mass distributions. 
We now compare the quantities extrapolated to the virial 
radius with predictions using both simulations and SDSS 

observations. 

Abundanc e matc hing studies like iMoster et al] (|2010h 
and iGuo et al] (|2010h make use of cosmological simulations 
and galaxy surveys to determine the mass dependence of 
galaxies and their preferred host haloes. The stellar mass 
enclosed within a 2R lclls aperture is known from our pop- 
ulation synthesis modelling, as shown in LFSF. The stellar 
mass profiles do not change significantly beyond 2_Ri en s . Thus 
we use the Af h aio-to-M ato i relation from [Moster et al] (|2010h 
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Sample 


Size 


Method 


M vir -range 


a 


C14 








[J.U M®\ 






B07 


39 


BCES 


0.06 - 20 


-0.199 ± 0.026 


9.12 ±0.43 


B07 


22 


BCES 


> 1 


-0.103 ±0.055 


7.71 ±0.58 


B07 


17 


BCES 


< 1 


-0.201 ±0.129 


9.46 ±2.11 


c 


18 


BCES 


0.004 - 0.4 


-0.401 ± 0.064 


7.03 ± 1.49 


c 


IS 


BS 


0.004 - 0.4 


-0.433 ± 0.078 


6.60+? 3 ,; 1 

—6.2 


c 


9 


BCES 


> 0.06 


-0.203 ±0.172 


16.98 ±12.87 


comb 


57 


BCES 


0.004 - 20 


-0.278 ±0.021 


9.62 ±0.41 


B07 


39 


BCES 


0.06 - 20 


-0.172 ±0.026 


9.0 ±0.4 




18 


BCES 


0.004 - 0.4 


-0.381 ±0.062 


12.02 ± 2.57 


CN7 


62 


N/K 


0.4- 100 


-0.15 ±0.13 


10.68 ±5.50 



Table 1. Slope a of the c — M vir relation with uncertainty for different samples, sample sizes, fitting methods and mass rang es (all errors 
la, ci 4 errors for BS-method are 68% conf. interterval around median). B07 denotes the sample of massive early types in lBuote et al.l 
(2007), c stands for the relation as in Fig. [4] comb gives the c ombin ed sample fit consisting of 39 objects of B07 and 18 lenses of the above 
cvir sample, CN7 are results from [Comerford &: Nat araian (20q3). The index to a sample name denotes concentrations normalized to 

2 = 0. 



to infer a virial mass and the corresponding scatter (taken 
at the la level). Note that the above abundance match- 
ing relations are based on Kroupa/Chabrier IMFs and thus 
consistent with the stellar masses used here. Fig. [5] shows 
the extrapolation of the NFW fits out to the virial radius. 
The figure also adds an extra point (r*",M*") marking 
the virial radius and mass from abundance matching. For 
10 out of the 18 lenses, M*, M turns out to lie within the 
la confidence region (grey shaded) around the original fit 
for AM(R). Accordingly, we show a further NFW fit (black 
curve) that is constrained to pass through (r*",M*"). 

In Fig. [6] we show how the c-M vir scatter plot changes 
when we impose abundance matching. The arrows point 
from values without the abundance-matching information 
to values including it, and the longer arrows are labelled. 
Going from Fig. [3] to Fig. [S] leads mostly to shifts along 
the direction of the relation, but the rms scatter with re- 
spect to the simple power-la w fit almost doubles , from 0.145 
to 0.258. In comparison, the lBuote et al.l (|2007l ) sample has 
a rms scatter of 0.180. A mildly increase d scatter can be 
found in simulations bv lShaw et al.l (|2006h for virial masses 
between ~ 3 x 10 13 and ~ 1O 15 M0, which is most likely due 
to the indistinguishability between substructure and main 
haloes. However, this cannot explain the increased scatter 
we find. We can conclude that the extrapolation to r vir in- 
ferred in Section [3] gives a reasonable extension to the c-M vir 
relation. Abundance matching, on the other hand, appears 
to introduce a large discrepancy in some cases. Three of 
the galaxies (MG2016, B1422 and B1030) show a shift to 
a much higher concentration when abundance matching is 
imposed. These are lenses for which M^f lies significantly 
below the extrapolation M vi[ . For MG2016, M* r M is even 
smaller than AM(< R) at the outermost radius of the lens 
model. Further three galaxies (Q0047, MG0414, SBS1520) 
show large shifts towards lower concentrations. They belong 
to the highest redshift galaxies in our sample and are probed 
in an exceptionally large radial range, up to 10% of the virial 
radius (see column 2i? lcns /r vir in Tab. [2]). Moreover, MG2016 
and SBS1520, which exhibit strongly discrepant M*", have 



reconstructed mass profiles with comparatively small uncer- 
tainties. 

So what is the reason for this discrepancy? A possi- 
ble explanation is suggested by a visible correlation between 
the length of the arrows and the environment of the lenses. 
For extrapolated virial masses much lower than M*" one 
may argue that lens profiles are shallower in group or clus- 
ter environments than in more isolated locations. This re- 
flects the inverse proportionality of concentration and en- 
closed mass and is a consequence of hierarchical structure 
formation. Extrapolating mass profiles from shallower pro- 
files leads necessarily to lower masses at r vir . In other words, 
if the M h aio(Af s toi) obtained from abundance matching is 
employed to determine Af*", we implicitly assume an iso- 
lated galaxy located within a "typical" halo with respect 
to its stellar content and the halo definition used in the 
abundance matching procedure. For lenses with extrapo- 
lated virial masses much larger than M*" the mere effect 
of the projected cluster environment might become more 
important, that is, although being relatively shallow, the 
projected total mass profile is strongly influenced by dark 
matter in the cluster acting as an additional convergence. 
This again causes the extrapolation to be significantly dif- 
ferent from M*". Examples for the latter case are MG2016 
and B1422, which are located in the densest environments 
among our lenses with large groups or clusters showing many 
nearby galaxies (cf. Table 1 in LFSF). 

All 8 lenses for which M*" is strongly discrepant with 
M vir are in dense environments, whereas for 6 out of 10 re- 
maining galaxies, there are no nearby objects reported so far. 
Current abundance-matching prescriptions do not consider 
environmental effects. Our results suggest that environment 
may significantly influence the M h aio(Af st oi) function. 



5 ADIABATIC CONTRACTION 

In the following section we will assess the extent to which 
the c-M vir r elation could be caused b y adiabatic contraction 
of the halo. iBlumenthal et all (| 1986T ) proposed that during 
the formation of galaxy-sized structures, the collapse of the 
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Figure 5. Extension of Fig. [3] to the virial radius and for the whole lens sample. In each panel there are three additional elements not 
present in Fig. \3\ the vertical dashed line marks the the r vir inferred from the NFW fit; the outer isolated point shows the virial radius 
and virial mass inferred from abundance matching; the black curve is an NFW fit constrained to go through the abundance-matching 
point. 



dissipative baryons towards the centre of the forming halo 
would extert a reaction on the dark matter density profile, 
making it steeper. This effect would mean that simple N- 
body simulations, such as those that led to the proposal of 



the NFW profile l|Navarro et alJll996f ). would underestimate 
the inner slope of the halo. 

The concentrations derived in Section would there- 
fore represent the state of the halo after adiabatic contrac- 
tion (AC). In this Section, we will refer to these concentra- 
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tions as Cf ina i. Before AC, the concentrations are thought to 
have a lower value c in i t . At present it is not clear wh ether 
or not Ci„it differed from c fina i (e.g.. IXbadi et aHl2010T ) and 
it is conceivable that the impact of AC on dark matter pro- 
files might be overestimated by commonly used recipes for 
baryonic cooling. Additional mechanisms suc h as dynamical 
encounters with smal ler structures (see e.g. |E1-Zant et al.l 
|2004| ; ICole et al]|201lf ) , or the e.jecta of b aryons triggered by 
supernovae- driven winds (|Larso nl ll974l ; iDekel fc Silkl 11989 ; 
iBrooks et al]|2007h could lead to the opposite effect, mak- 
ing the inner slope of the dark matter profile less cuspy. In 
this paper we only consider the effect from the more funda- 
mental process of contraction during the formation of the 
halo. 

5.1 Comparing Adiabatic Contraction Models 

To analyze this issue, we make use of the halo contrac- 
tion program of iGnedin et al.l (|200J), which computes the 
change in the dark-matter density profile under AC, keep- 
ing rM(< r) conserved. To take account of a wide range of 
orbit eccentricities the code invokes the power-law 

f/r vir = A(r/r VII ) w (12) 

to describe the mean relation between orbit-averaged and 
current radius, and modifies the adiabatic invariant to rM(< 
f). Eg. 1121 changes the eccentricity distribution of the mass 
profile, which is thus distorted by the usage of a mean radius 
in the invariant. Parameter A defines the maximum eccen- 
tricity and causes rM(< f) to be larger than rM{< r) for 
r/r V j r < 0.44 and smaller for r/r vir > 0.44. A larger invari- 
ant means more mass in the center at the expense of the 
outer parts of the halo. The parameter w defines how strong 
the shift is. The smaller w the fewer mass at the center. 

The case A = w = 1 therefore co r respo nds to the orig- 
inal prescription of iBlumenthal et al.l fl986,), where the or- 
bits are assumed to be completely circular. This case can 
be understood as an upper limit to AC. The program pro- 
vides the necessary resolution for comparison with our data, 
i.e., down to 10 _3 r vir . The input parameters are the 
baryon fraction enclosed within r vir , the baryon scale length 
and the initial concentration of the dark matter halo, c init . 
We take the baryon fraction as M ste i(< 2i? lcns )/M vir , where 
M Bte i(< 27? lcns ) denotes the stellar mass enclosed in the total 
reconstructed radial range. For the baryon scale length we 
use the fitted Hernquist scale radius derived in Section [3] 
This is preferred to making use of R a w th X 1.8 (|Hernquistl 
1990), because our measured R a - derived from the Petrosian 
radius - do not agree precisely with the half-light radius of 
Hernquist profiles, which is a consequence of projected radii 
and circularized mass profiles. Furthermore, the Hernquist 
profile is originally used for the surface brightness distribu- 
tion, whereas we fit in this case surface mass profiles. The 
r^/iZc best-fit values turn out to be mostly lower but close 
to 1/1.8. 

We run the contraction routine for a grid of pa- 
rameters ( Ci n it, fb, r h /r wil ) ranging from (5,0.005,0.001) to 
(60,0.135,0.015) in steps of (1,0.01,0.002). We then fitted 
the contracted profiles, via Eq. [5] to the data AM(< R) for 
R/r v -„ ranging from ~ 0.006 to 0.12. 

There are a number of uncertainties entering the anal- 
ysis: 




0.00 0.03 0.06 0.09 0.12 

fb 



Figure 8. Final concentration versus baryon fract ion depending 
on siz e of radial window for the AC prescription of Gned in et al.l 
(2004). Bright to dark grey corresponds to aperture sizes (~ 
0.003 - 0.06)r vir ,(~ 0.005 - 0.09)r vir and (~ 0.006 - 0.12)r vir . 

(i) since the radial extent of a reconstructed profile is lim- 
ited to 2 times the angular i?i on5 and a finite resolution, the 
aperture size changes from lens to lens, 

(ii) in order to mimic the limited probed range (hence- 
forth called aperture) by an equivalent range in the con- 
tracted profile, -Rions must be expressed in units of r viT , 

(iii) baryon fraction as well as baryonic scale length de- 
pend on M vir and r vir which are extrapolated quantities with 
their respective uncertainties. 

Figure [7] shows the results for three scenarios of adia- 
batic contraction. The leftmost panels show initial versus 
final (i.e. contracted) halo concentrati on. The top row 
corres ponds to the original proposal of IBlumenthal et all 
(1986) (A = 1.0, w = 1.0, i.e., no correction for anisotropic 
orbits). In this case, we illustrate the increase in concentra- 
tion as a blue line for fixed values of fb and rd/r vir . The 
growing concentration towards low Ci nl t is a consequence 
of the interplay between radial aperture - i.e. the extent 
of the extrapolation - and the region affected by AC. 
The smaller the initial concentration the larger CH na i/cinit 
towards small radii for the same fb (we refer to this as 
the AC-sensitive case.). As Ci nl t increases the difference 
between final and initial profile subsides. Note that in 
our analysis, the further out we can probe the halo, the 
less affected is the fit and the extrapolation. However, for 
different combinations of fb and rh/r vir , similar curves fill 
the grey-shaded region. To enable proper differentiation 
with respect to initial concentration, we choose different 
shades of gray. The middle (rightmost) panels show the 
final concentration versus baryon fraction (baryon-to-virial 
radius fraction). The gray shaded regions map the same 
areas as those in the leftmost panels. T he black dots with 
error bars represent our data. For the IBlumenthal et all 
(|l986h case (top), there is clearly a disagreement between 
observationally inferred and contracted profiles. Especially 
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0.005 0.010 0.015 



Figure 7. Adiabatic contraction following IGnedin et alJ (2004). Each row corresponds to a different set of values for the parameters A 
and w in Equation 112> . Each row shows the final concentration against: (left) the initial concentration, (middle) the baryon fraction, 
and (right) the scale radius of the baryons. The initial concentration can be read off the middle and right panels by comparing the gray 
bands with the left panel of the same row. The colored lines illustrate the effect of holding fixed two out of (i) initial concentration, 
(ii) baryon faction and (iii) baryon scale radius. 



the low- ft and low-rj,/r vir regions show significant de- 
parture from even the lowest final concentrations of the 
generic haloes. The middle row of F igure [7] shows the AC 
prescription of IGnedin et al.l (|2004T ). that implements eq. 
1121 with fiducial values A — 0.85 and w — 0.8 to take 
into account eccentric orbits. This phenomenologically 
motivated ansatz leads to smaller concentrations. There is 
still significant disagreement between data and simulated 
contraction. The behavior of Cvirjinal as a function of 
rh/r V ir for constant c vir j ina i and ft is indicated by the 
solid magenta line. For the panels in the bottom row of 
Fig. we changed the pre-defined values of A and w to 
0.42 and 0.4, respectively. The orange line shows for fixed 
c init and rh/r vir the final concentration as a function of /;,. 



These values for A and w give good agreement with the 
lensing data even for low c init , between 1 and 10. Compared 
to the AC prescriptions shown in the top and middle rows, 
the range of final concentrations is narrower, corresponding 
to a shallower Cnnia-fb relation (middle panels). The latter 
can be equivalently expressed in terms of mass not drawn 
into the central region < 0.1r vir . Comp ar ing m ass profiles 
contracted according to IGnedin et al.l (|2004l ) with the 
(A — 0.42, w — 0.4)-case shows that the latter transports 
less mass (~ 0.4% of the virial mass) into the central halo 
region, < 0.1r vir . 



One of the intriguing results of this study is that even 
with a simple assumption of a common, mass-independent 
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initial concentration, most of the final concentrations can 
be explained if (A,w) are conveniently adjusted and th is 
allowed to vary within uncertainties. There are a variety of 
results we summarize in the following. 

• There is slight evidence for lenses with lower baryon 
fraction to require higher initial concentrations. 

• Both smaller c ini t and smaller rh/r viT produce steeper 
Cfi„ai(/t) curves. This effect is independent of the AC- 
sensitive case at very low c init (explained above). 

• When A and w are reduced, ea na i(/f)) and Ca n!1 \(rh) be- 
come flatter, i.e. the differently shaded c init regions of the 
leftmost panels are mapped to narrower regions in the mid- 
dle and right panels. Moreover their overlap is reduced. 

To study how sensitive our results are to the radial 
range, we additionally provide in Fig. [5] the results for 
reduced aperture sizes, i.e. (~ 0.005 — 0.09)r vir and (~ 
0.003 — 0.06)r vir using the parameters A = 0.85 and w = 0.8. 
We find that even for the smallest radial range fits do not 
yield agreement with our lens data. The case that we un- 
derestimate the radial extend of our lenses by a factor of 
two is in light of the relatively small uncertainties of the 
mass profiles and the errors attached to r s unlikely. Larger 
radial apertures yield even more disfavored final concentra- 
tions compared to our lens data. 
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Figure 9. As in Fig. [4] but with initial concentrations. The col- 
ors of the columns from dark to bright correspond to the 68%, 
95% and 99% confidence intervals of a range of (c,/t,,r h /r vir ) val- 
ues which produce a Cf ina i in agreement with our data. For com- 
parison we include Cf ina i(M vir ) and results from simulations by 
IMaccio et all (I2007T) . The solid line indicates their mean concen- 
tration together with a 2a band (grey region). 



5.2 Initial Concentration from Weak Adiabatic 
Contraction 

For the weak AC case (A = 0.42, w — 0.4) we compare 
the number of different (c iait ,rh/r vil ) combinations produc- 
ing final concentrations in agreement with our lens data and 
infer a Ci„it-M vir plot as before, enriched by the information 
of the frequency distribution of initial concentrations (see 
Fig. [9]). Although most of the data can be reproduced even 
by few initial concentrations of ~ f , most of the (c ia i t ,rh/r vil ) 
combinations with high c init produce final concentration in 
agreement with c fina i and /;,. The hue of the magenta column 
indicates the frequency distribution of c init values whereas 
the 68% (99%) confidence interval is highlighted by strongest 
(faintest) color. 

Certainly, no strong quantitative conclusions can yet be 
drawn from these results, but judging by the confidence re- 
gions a strongly flattened c in it(M vir ) relation seems likely. 
These results are in agreement with the i nital halo concen- 
trations estimated bv lLintott et al.l (|2006T l. who used a sim- 
ple model of spherical collapse. In that model, massive galax- 
ies from density fluctuations between 2 and 3<r - roughly 
mapping the same mass range as our lensing galaxies - were 
found with initial concentrations between 3 and f 0, with the 
most massive ones having lower concentrations. A flattened 
low-mass end of the c — M vir relation is expe cted by simu- 
lation s. To show this we include results from IMaccio et al.l 
(2007) (solid line and gre y 2a band in Fig. c (M, z) curves 
based on a toy model bv lBullock et al.l l|2001f ) for redshifts 
in a range from to f .4 are in good agreement with results 
from N-body simulations. The toy model includes the free 
parameter K which takes into account the contraction of the 
inner halo beyond that required by the top-hat formation 
scenario. This contraction parameter is fixed for all haloes in 
their simulation. The difference between the simulated and 
observed c — Af vir relation is a well-known issue and matter 



of ongoing studies. It is however worthwhile to mention that 
this discrepancy is even stronger for low virial mass. From 
the comparison between c ini t-M vir found in this study and 
simulations that investigate the redshift-dependence, we can 
conclude that Adiabatic Contraction alone is not enough to 
explain the slope of the relation. 



6 CONCLUSIONS 

Strong gravitational lensing on galaxy scales constitutes a 
powerful tool to characterize dark matter haloes. In addi- 
tion, combining photometric studies with stellar population 
synthesis allows us to assess the interplay between baryons 
and dark matter in the central regions of galaxies. This pa- 
per extends the work of LFSF by exploring in detail the 
concentration of the dark matter haloes of 18 massive early- 
type lensing galaxies. On a concentration-virial mass dia- 
gram (Fig. 2| we find these haloes to confirm and extend 
towards lower mas ses the relationship observed in X-rays by 
iBuote et al.l (|2007T l. 

Our sample includes information about the baryon frac- 
tion, enabling us to explore the validity of adiabatic con- 
tractio n p rescriptions, such as the one of iBlumenthal et al.l 
l|l986h or iGnedin et all (|2004T ). We find that the standard 
modelling gives rather high final concentrations compared 
to our observations ( Fig. O. A tweak of t he parameters in 
the AC prescription of lGnedin et al.l (|2004T ) cause the gain in 
mass of the central region (< 0.1r vir ) to be ~ 4 x 10" 3 Af vir 
lower than in the previous case, which helps to solve the 
discrepancy. Furthermore, this results in a rather flat re- 
lationship between initial concentration (i.e. pre-AC) and 
halo. 

We emphasize that this paper focusses on the funda- 
mental aspect of adiabatic contraction caused by the col- 
lapse of the baryons during the formation of the haloes. 
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Additional m echanisms acting later, arising from d ynamical 
interactions (|E1-Zant et all |2004| ; ICole et all l201ll i or stel- 
lar feedback resulting in the expulsion of baryons (see e.g. 
iRead fc Gilmorell2005l ; iBrooks et al.ll2007l ) may alter the in- 
ner slope of the dark matter halo, although these mech- 
anisms are expected to be more important in lower mass 
galaxies. Nev e rthele ss, the tweak in the AC prescription of 
iGnedin et al.l |2004!) could be interpreted as one of these 
mechanisms playing a role in the evolution of the structure 
of the haloes. In any case, our analysis suggests that adia- 
batic contraction can explain only part of the c-M vir trend 
but is unlikely to be the sole origin of it. 
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12 D. Leier, I. Ferreras, P. Saha 



Lens z L r s p s r vir - r h ^ M vir 
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Table 2. Lens, redshift of the lens zl, NFW scale radius r s , NFW scale density p B , inferred virial radius r vir , outermost radius of the 
mass profiles in terms of virial radius 2i?i ons /r vir (for the innermost radius multiply by 1/19), Hernquist scale radius r^, Hernquist scale 
luminosity M st ei/Y, the virial mass M vir as defined in Equation 151 the concentration as defined in Equation [9] 
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